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ABSTRACT 



We study two samples of local galaxies, one is UV (GALEX) selected and the 



other FIR (IRAS) selected, to address the question whether UV and FIR surveys 

see the two sides ('bright' and 'dark') of the star formation of the same population 

of galaxies or two different populations of star forming galaxies. No significant 
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difference between the L tot (= L 60 + Lpuv) luminosity functions of the UV and 
FIR samples is found. Also, after the correction for the 'Malmquist bias' (bias for 
flux limited samples), the FIR-to-UV ratio v.s. L tot relations of the two samples 
are consistent with each other. In the range of 9 < \og(L tot /L Q ) < 12, both can be 
approximated by a simple linear relation of \og{L^/ Lpuv) = log(L toi /L ) — 9.66. 
These are consistent with the hypothesis that the two samples represent the same 
population of star forming galaxies, and their well documented differences in L tot 
and in FIR-to-UV ratio are due only to the selection effect. A comparison between 
the UV luminosity functions shows marginal evidence for a population of faint 
UV galaxies missing in the FIR selected sample. The contribution from these 
'FIR-quiet' galaxies to the overall UV population is insignificant, given that the 
K-band luminosity functions (i.e. the stellar mass functions) of the two samples 
do not show any significant difference. 

Subject headings: dust: extinction - galaxies: luminosity function, mass function 
- infrared: galaxies - ultraviolet: galaxies 



1. Introduction 

The evolution of star forming galaxies tells much about the history of the universe. 
The star formation activity in these galaxies can be best studied by observing the emission 
from young massive stars in the rest frame UV and FIR. The UV observations record the 
direct light from the hot young stars, and the FIR observations collect star light absorbed 
and then re-emitted by the ubiquitous dust. A complete picture of star formation in the 
universe can only be obtained when the observations in these two wavebands are properly 
synthesized. Indeed, our knowledge on the star formation history of the universe has been 
mostly derived from deep surveys in the rest frame UV and FIR. Many studies have been 
devoted to methods of deriving star formation rate of individual galaxies using the UV or 
FIR luminosities (Calzetti 1997; Meurer et al. 1999; Buat & Xu 1996; Iglesias-Paramo et al. 
2005), and the strengths and shortcoming of these methods have been discussed thoroughly 
in the literature (Kennicutt 1998; Adelberger & Steidel 2000; Bell 2002; Bell 2003; Buat 
et al. 2005; Kong et al. 2004; Iglesias-Paramo et al. 2006). However, an arguably more 
important issue is the selection effect of the surveys that can be summed up by the following 
question: Do UV and FIR surveys see the two sides ('bright' and 'dark') of the star formation 
of the same population of galaxies, or do they see two different populations of star forming 
galaxies? This is important because if the correct answer is the latter, then even if one 
can estimate accurately the star formation rate for galaxies in surveys in one band, the star 



-3- 



formation in galaxies detected in the other band is still missing. Actually this question is in 
the core of an on-going debate on whether the SFR of z~ 3 universe can be derived from 
observations of Lyman-break galaxies, which are UV selected star forming galaxies at z~ 3 
(Adelberger and Steidel 2000), given that SCUBA surveys in sub-millimeter (rest frame FIR 
for z> 2) detected many violent star forming galaxies at about the same redshift that are 
not seen by LBG surveys (Smail et al. 2001; Smail et al. 2004). 

There have been limited overlaps between rest frame UV surveys and rest frame IR 
surveys. In the SCUBA survey of LBG's (Chapman et al. 2000), only one LBG was detected. 
As summarized in Adelberger and Steidel (2000), only a couple of SCUBA sources are bright 
enough in optical to be detected in LBG surveys. The situation is better for z~ 1 star 
forming galaxies, which now can be routinely identified by large scale spectroscopic surveys 
and multi-band optical surveys, and which have been detected in abundance in mid-infrared 
by ISOCAM deep surveys (Elbaz et al. 2002; Hammer et al. 2005) and Spitzer surveys 
(Le Floch 2005), and in UV by GALEX (Arnouts et al. 2005; Burgarella et al. 2006). 
However the extrapolation from the mid-IR to the total dust emission is very uncertain and 
may subject to significant evolution itself. The same criticism can also be applied to the 
comparisons between rest frame UV and MIR sources at z~ 2, the latter being detected 
recently by Spizter at 24/mi (Chary et al. 2004; Shupe et al. 2005). Because of the relatively 
high confusion limits for surveys in the Spitzer MIPS 70fim and 160/im bands, thorough 
comparisons of rest frame UV and FIR sources of z> 1, down to luminosity levels fainter 
than the 'knee' of the luminosity functions of both bands, may have to wait until the launch 
of Herschel (Pilbratt 2005). 

In this paper, we investigate the difference and overlaps of the UV and FIR selected 
samples in the local universe in an attempt to shed light on the selection effects of high-z 
samples similarly selected. The UV data are taken from observations by Galaxy Evolution 
Explorer (GALEX) and the FIR data are taken from the IRAS database. Several papers have 
been published using these data. Martin et al. (2005, hereafter M05) derived the local (z=0) 
bivariate luminosity function for the FUV (1530A) and FIR (60fim) bands, which shows that 
the FUV luminosity saturates at about 2 10 10 L Q while the FIR luminosity can be as high as 
~ 10 13 L Q . This is consistent with a very strong dependence of the FIR/FUV ratio on the 
total luminosity (L tot = L F i R + L FU v)- The luminosity function of L tot has a log-normal form. 
Buat et al. (2005) compared the extinction properties of local UV and FIR selected galaxies 
and found that the mean NUV (2267A) extinction of UV selected galaxies is significantly 
lower than that of FIR selected galaxies (~ 1 mag vs. ~ 2.5 mag). Iglesias-Paramo et al. 
(2006) carried out an extensive study, using combined GALEX and IRAS data, on the UV 
and FIR emission as star formation indicators and found a rather modest star formation 
activity for local star forming galaxies. Pre-GALEX studies on comparisons between UV 
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and FIR selected samples can be found in Buat & Burgarella (1998), Buat et al. (1999) 
and Iglesias-Paramo et al. (2004). In this work, we will study some statistics free of the 
selection effect in order to check quantitatively how much the UV and FIR selected samples 
differ/overlap with each other. The paper is arranged as following: After this introduction, 
the data sets analyzed in this paper are presented in Section 2. Major results are listed 
in Section 3. Section 4 is devoted to the discussion. Through out this paper, we assume 
tt A = 0.7, Vt m = 0.3, and H = 70 km sec" 1 Mpc" 1 . 



2. Data 

The data sets are basically the same as those in Iglesias-Paramo et al. (2006) and Buat 
et al. (2005). The original UV selected sample (Iglesias-Paramo et al. 2006) includes 95 
galaxies brighter than NUV = 16 mag selected from GALEX Gl stage All-sky Imaging Survey 
(AIS), covering 654 deg 2 . From these we exclude 1 galaxy, 2MASX_J20341333-0405, which 
does not have measured redshift. The FIR selected sample is also taken from Iglesias-Paramo 
et al. (2006). From the original sample, including 163 galaxies with / 60 > 0.6 Jy in 509 deg 2 
sky covered both by GALEX AIS and IRAS PSCz (Saunders et al. 2000), 2 are excluded: 
NGC 7725 (no redshift) and IRAS-F00443+1038 (not a galaxy). Consequently, the UV and 
FIR selected samples studied in this paper have 94 and 161 galaxies, respectively. The K s 
(2.16//m) band magnitudes K to t are taken from Extended Source Catalog (XSC) of 2MASS 
(Jarrett et al. 2000). Since this is very close to the classical K (2.2/zm) magnitude, we will 
call it K magnitude hereafter for the sake of simplicity. For both the UV selected sample 
(94 galaxies) and the FIR selected sample (161 galaxies), each has 12 galaxies undetected 
by 2MASS. According to the sensitivity limit of 2MASS XSC (Jarrett et al. 2000), upper 
limits of K=13.5 mag are assigned to the undetections. 

Morphological classifications were searched in NED. For those galaxies without mor- 
phological classification in the literature, images taken from SDSS, DSS and 2MASS (in 
the order of priority) were inspected and eye-ball classification was carried out. Galaxies 
included in Atlas of Peculiar Galaxies (Arp 1966), Southern Peculiar Galaxies and Associa- 
tions (Arp & Madore 1987), and Catalog of Isolated Pairs of Galaxies (Karachentsev 1972) 
are classified as peculiar, interacting or mergers. For a few galaxies that are faint (b > 15 
mag) and small (< 10") the classification can be very uncertain. Most of such galaxies are 
in the FIR selected sample, and very often there is clear sign of interaction (close companion 
of similar brightness and/or diffuse tidal features). In Fig.l the distributions of morpho- 
logical types (not including QSOs and ellipticals) of the two samples are compared. The 
overall overlap between the two distributions is about 60%. There is a significant excess 
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Fig. 1. — Morphological type distributions of the UV and FIR selected samples. 
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of Pec/Int/merg galaxies in the FIR selected sample (39%) compared to those in the UV 
selected sample (14%). For normal galaxies both UV and FIR selected samples peak in the 
bin of Sab/Sb/Sbc. Detailed analysis shows that the median type for normal UV galaxies 
is Sc and that of normal FIR galaxies is Sb. The FIR selected sample is tilted toward the 
earlier spirals whereas the UV sample has more late type (later than Sc) galaxies. 



3. Results 

3.1. Comparisons of luminosity functions of UV and FIR galaxies 

Much of the difference between UV and FIR selected samples can be traced back to 
a single selection effect: UV observations detect preferentially galaxies with low L FIR /L UV 
ratios, and in contrast FIR observations select galaxies with high Lfir/Luv ratios. Since 
the FIR/UV ratio is a good indicator of dust attenuation (Xu & Buat 1995; Buat & Xu 
1996; Meurer et al. 1999; Gordon 2000), it follows that UV samples select galaxies with 
significantly lower dust attenuation than galaxies in the FIR selected sample: Buat et al. 
(2005) found a median FUV attenuation of A(FUV) = 0.81^3 mag for the UV selected 
sample, compared to a A(FUV) = 2.llJ'g mag for the FIR selected sample. 

It has been well established that there is a strong correlation between luminosity and 
dust attenuation in the sense that more luminous galaxies have higher dust attenuation 
(Wang & Heckman 1996; Buat & Burgarella 1998; Adelberger & Steidel 2000; M05). Fig.2 
shows that galaxies in both FIR and UV samples follow the strong L 60 /Luv vs. L tot cor- 
relation. On the other hand, UV galaxies in general have significantly lower L toi and lower 
L 60 /Luv ratios for a given L tot compared to FIR galaxies (Buat & Burgarella 1998; Adel- 
berger & Steidel 2000; Iglesias-Paramo et al. 2006). Can these trends be attributed solely to 
the selection effects, or do they reflect some intrinsic differences between the two populations? 

In order to answer this question, we have to compare the statistics of L tot and of Lgo/Luv 
that are free from the selection effect. The selection effect is introduced by the so-called 
'Malmquist bias' on both flux limited samples: For a given L ioi , galaxies with higher FIR- 
to-UV ratios have brighter Lqq, therefore can be seen at larger distances (i.e. having a larger 
maximum finding volume V max ) in a f^ limited sample. Similarly, for a given L tot , galaxies 
with lower FIR-to-UV ratios have higher Lyy and therefore larger V max in a UV flux limited 
sample. In what follows we shall compare the L tot luminosity functions (LFs hereafter) of 
the two samples to exam whether they have the same intrinsic L tot distributions. Because 
LFs are luminosity distributions of galaxies in a unit volume, they are not subject to the 
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Fig. 2. — The Lq /L fuv ratio vs. L tot (Lq +L fuv ) plot for UV (blue symbols) and FIR (red 
symbols) selected galaxies. The cosmic mean of the FIR/UV ratio, < Pdust/ Pfuv > is taken 
from Takeuchi et al. (2005), assuming L 60 = 0.4 x L^ust- 



bias discussed above. 

Here we exclude the sources whose IRAS fluxes are affected by the cirrus. Also UV 
galaxies not covered by the IRAS survey are dropped. This reduces the FIR sample to 
151 galaxies and the UV sample to 81 galaxies. The IRAS detections of 5 UV galaxies 
are confused with other UV sources, therefore the corresponding IRAS fluxes are treated 
as upperlimits. Altogether 14 UV galaxies have only upperlimits for the IRAS flux. For 
galaxies in the FIR selected sample, 14 have only upperlimits for the FUV flux. 

Define 0f o Y v (L k ) as the L tot LF of UV selected galaxies at log(L to t) = Lk, 0fuv(Li) 
the FUV (1530A) LF at log(Lpuv) — Li, and P kj j the conditional probability of finding UV 
galaxies of log(L FU v) = L ; in the bin of L k — 0.55 k < log(L tot ) < L k + 0.55 k . Then 



/.FUV 
y tot 



(L k ) = ^P kii FUV (L i )<5 i /5 k . (1) 



Similarly, the L tot luminosity function of FIR selected sample can be derived using the 
formula: 

OW) = J>k'j<MLj)V^ (2) 

j 
where P k ',j the conditional probability of finding FIR galaxies of logLgo = Lj in the bin of 
L^ — 0.55y < log(L tot ) < L k ' + 0.55k'. Data in our two samples are used in the calculations 
of the conditional probability functions P k ,i and P k ',j- In order to take into account the 
information content in the upper limits, the Kaplan-Meier (KM) estimator (Kaplan & Meier 
1958; Feigelson & Nelson 1985; Schmitt 1985) has been applied in these calculations. We 
have chosen 5{=1 mag for the Lpuv bin width, Sj— 0.5 dex for the L@o bin width, and 5^=0.5 
dex for the L tot bin width. Other choices of the bin widths result in LFs with either larger 
scatters (bin widths too narrow) or coarse resolutions (bin widths too broad). The FUV LF 
and Leo LF are taken from Wyder et al. (2005) and Takeuchi et al. (2003), respectively. 

The results are listed in Table 1 and plotted in Fig.3. In the L tot range where they 
overlap, the LFs of the two populations are consistent with each other. The solid line is 
the best fit of the L tot LF of M05, derived from a combined sample of UV and FIR selected 
galaxies. It is a log-normal function with the center at log(L toi /L ) = 7.43 and o = 0.87. 
In bins of log(L tot /L ) > 10, our LFs are marginally higher than that of M05. In order to 
check whether this indicates over-estimation in our results, we also compared with the Lqq 
LF of Takeuchi et al. (2003). There is a good agreement between our results and that of 
Takeuchi et al. (2003) for bins of log(L tot /L Q ) > 11 (where L 60 always dominates L tot ), both 
are slightly higher than that of M05. At log(L toi /L ) = 9, our results for both samples are 
lower than that of M05, possibly due to uncertainties caused by the small size of our samples 
compared to that of M05. 
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Fig. 3. — The L tot (L 60 + Lpuv) luminosity functions of UV galaxies (open diamonds) and 
FIR galaxies (solid squares). 
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The above result is consistent with that the UV and the FIR samples are drown from 
the same population of star-forming galaxies. However, the L tot LF comparison could be 
insensitive to some differences. For example, in bins where L tot is dominated by L 60 , the dif- 
ferences between ~Lfuv distributions of two samples can be hidden by the similarity between 
Leo distributions, and vise versa. Therefore, in what follows we shall calculate the Leo LF 
of UV galaxies and compare it with the L@o LF of IRAS galaxies (Takeuchi et al. 2003), 
and calculate the Ijfuv (1530 A) LF of FIR galaxies and compare it with that of GALEX 
galaxies (Wyder et al. 2005). 

The formalism for the calculations of Lgo LF of UV galaxies, (^^(Lqo), and of ~Lfuv 
LF of FIR galaxies, 0puv(Lfuv), is the same as that used in the calculations of L tot LFs. 
One only needs to replace L tot by L 60 in Eq(l), and by L FUV in Eq(2). 

The results are listed in Table 2 and Table 3. In Fig. 4, <f)Q^ v (Lj) is compared with 
the 60 fim luminosity function of IRAS sources (Takeuchi et al. 2003). It appears that UV 
galaxies can account for the FIR luminosity function up-to Leo ~ 10 11 ' 5 L Q . Only ULIRGs 
of L 60 ^ 10 12 L Q are missing in the UV sample. This is because ULIRGs are very rare in the 
local universe, and they are much fainter in UV. Therefore they are probed by UV surveys 
in a very much smaller volume compared to that probed by the FIR surveys. It should be 
pointed out that the UV LF of Wyder et al. (2005) excludes the contribution from broad-line 
AGNs identified using SDSS spectra. These are UV/optical selected QSOs and Seyfert 1 
galaxies. According to Sanders et al. (1989) and Spinoglio & Malkan (1989), these sources 
never contribute more than 10% of the IR LF in the whole range of FIR luminosity. The 
comparison between FUV luminosity function of the FIR selected sample and the GALEX 
FUV luminosity function (Wyder et al. 2005) is in Fig. 5. It shows that UV galaxies brighter 
than L*(FUV) (~ 10 9 ' 5 L ) are fully represented in the FIR selected sample. In fact there 
is a significant excess in the brightest bin (Mpuv — ~~ 21) of the UV LF of FIR sources 
compared to the UV LF of Wyder et al. (2005), likely being caused by the exclusion of 
the broad-line AGNs in the latter. There is a marginal evidence for fainter UV galaxies 
of Lpuv < 10 9 ' 5 L Q being under-represented in the FIR selected sample, suggesting that a 
population of 'FIR-quiet' UV galaxies might be missing in the FIR selected sample. 



3.2. FIR-to-UV v.s. L tot relations of UV and FIR galaxies 

Let R = log(L 60 /L F uv) = log(L 60 ) — log(L F uv)- For UV galaxies with a given L to t = L k , 
a 'Malmquist-bias-free' (i.e. selection effect free) indicator of mean FIR-to-UV ratio can be 
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Table 1. The L toi (L@o + Lfuv) Luminosity Functions of UV and FIR Selected Galaxies. 
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Table 2. The FIR (60//m) Luminosity Function of UV Selected Galaxies. 
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Fig. 4. — The Lfir (60/im) luminosity function of UV selected galaxies compared to the 
IRAS 60 fim luminosity function (Takeuchi et al. 2003). 
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Fig. 5. — The FUV (1530A) luminosity function of FIR selected galaxies compared to the 
GALEX FUV luminosity function (Wyder et al. 2005). 
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defined as follows: 

where Py is the conditional probability of finding UV galaxies of log(LFuv) — Li i n the FIR 
luminosity bin Lj — 0.56j < log(L 60 ) < Lj + 0.5<5j, and the summation goes through both 
indexes i and j including all bins satisfying the condition Lk — 0.55k < log(10 Li + 10 Lj ) < 
Lk + 0.55k- A similar FIR-to-UV ratio indicator can be defined for FIR galaxies: 

MU) " E^MW ( 4) 

where Py is the conditional probability of finding FIR galaxies of log(Lgo) = Lj in the FUV 
luminosity bin L; — 0.55; < log(L FUV ) < L ; + 0.55;. The variance of Ruv(Lk) and that of 
R-Fm(Lk), respectively, are: 

^Uv( L k) = — ^ p T 7TTT , (5) 



and 



2 , T , Ej.iKLj " Li) - RmMPPy^oMij ,^ 

<W L k) = — ^ , , u • (6) 



Results for Ruv, Rfir, o"uv, and ctfir, are listed in Table 4. As shown in Fig.6, there 
is no significant difference between the Ruv v - s - L tot relation of UV galaxies and the Rfir 
v.s. L tot relation of FIR galaxies, again in consistence with the hypothesis that the two 
samples represent the same population, and their difference in Fig2. is due to the selection 
effect. Both R v.s. L ioi relations can be approximated by a simple linear relation: R = 
log(L tot ) — 9.66, as shown by solid line in Fig.6. In the L tot range covered by our samples, 
this relation is slightly lower than the non-linear relation between R and log(L tot ) (dashed 
curve in Fig.6) derived by M05 from the bi-variate function of their combined sample. It 
should be pointed out that the simple linear relation should not be extrapolated to galaxies 
of L to t ;$ 10 9 L Q , where a flatter relation is more likely (M05). 



3.3. K band luminosity functions and stellar mass distributions 

The NIR K band luminosity, very insensitive to both the dust extinction and the star 
formation history variation (Bell & De Jong 1991; Bell et al. 2003), is the best stellar 
mass indicator. The stellar mass distribution is one of the most important characteristics 
defining galaxy populations, therefore we would compare the K band LF of UV galaxies 
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Table 3. The FUV (1530 A) Luminosity Function of FIR Selected Galaxies. 
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Table 4. The FIR-to-UV ratio (R = log{L m /L FUV )) v.s. L tot (L 60 

UV and FIR Selected Galaxies. 
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Fig. 6. — The FIR-to-UV ratio (R = log(L 60 / Lpuv)) v - s - L iot (L 60 + Lpuv) relations for UV 
and FIR Selected Galaxies. Solid line: R = log(L tot ) — 9.66. Dashed curve: R v.s. log(L tot ) 
relation derived by Martin et al (2005) for a combined sample. 
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Fig. 7. — K band luminosity functions (stellar mass distributions) of UV and FIR selected 
samples. 
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with that of FIR galaxies. Because of the present of upper limits in the K band fluxes in 
both UV and FIR samples, we exploit the same formalism as presented in Eq(l) and Eq(2), 
and use KM estimator in calculating the conditional probability functions P(M — 0.55 < 
M K < M + 0.5<5|L FUV ) and P(M - 0.55 < M K < M + 0.55|L 60 ). In Fig. 7, the resulted K 
LFs of the two samples are compared with each other. No significant difference is found 
between them. It is interesting to note that both K LFs are consistent with the K LF 
of late-type galaxies derived by Kochanek et al (2001), specified by a Schechter function 
with 0o = 0.0101, a = -0.87 and M, = -22.98 + 5. * alogl0(h ) - 5, where h = 0.7 and 
5 = 0.2 (the difference between the isophotal magnitude and the 'total' magnitude, Cole 
et al. 2001). The conversion factor M stars /LK = 1.32M Q /L , which is derived for a stellar 
population with constant star formation rate and a Salpeter IMF (Cole et al. 2001), is 
assumed when converting the K band luminosity to stellar mass. 



4. Discussion 

Our results indicate that bulk of z=0 galaxies selected in the UV and FIR samples 
are from the same population of active star forming galaxies. In particular, galaxies in the 
two samples have indistinguishable L tot LFs. And their FIR-to-UV ratio v.s. L tot relations, 
after correction for the Malmquist bias, are consistent with each other. Therefore, the well 
documented results that galaxies in the UV flux limited samples tend to have lower L tot and 
lower FIR-to-UV ratios for a given L toi than those galaxies in the FIR flux limited samples 
are purely due to the selection effect. 

The only sign for a possible difference between UV and FIR populations is a marginal 
deficiency of galaxies of low UV luminosity in the FIR selected sample, indicating the exis- 
tence of an 'FIR-quiet' UV population. Indeed, it has been known that there is a population 
of low-metallicity, low dust content 'blue compact dwarf galaxies. The prototype is I ZW 18, 
the galaxy with one of the lowest metallicity of 1/50 solar (Searle & Sargent 1972). I ZW 18 
has never been detected in FIR. The FUV magnitude of I Zw 18 derived from its GALEX 
image is 15.75 mag (Gil De Paz, private communication). Its IRAS upper limit of feo^m < 0.2 
Jy corresponds to a upperlimit of L 60 /Lpuv < 0.27. There are only a few percent of galaxies 
in our UV sample have such low Lqo/Lfuv ratio, indicating a low contribution from these 
'FIR-quiet' galaxies to the overall UV population. This is in agreement with the result in 
Fig. 7 which shows no significant difference in the K LFs (i.e. stellar mass functions) of the 
UV and FIR galaxies. Furthermore, because they have rather low UV and FIR luminosities, 
these galaxies contribute negligibly to the total SFR of the local universe. It will be inter- 
esting to know whether in the earlier universe more star forming galaxies becoming 'FIR 
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quiet', given the lower metallicity in high z galaxies and marginal evidence for net increase 
of the faint end of the UV LF (Arnouts et al. 2005). The new results of Burgarella et al. 
(2006) on LBG galaxies at z~ 1 suggest the existence of a population of low-attenuation, 
bright UV galaxies at that redshift. 

There are no ULIRGs in our UV sample. It is generally true that ULIRGs are absent in 
UV samples of sizes less than a few 1000s. In the local universe, LIRGS/ULIRGs contribute 
less than a few percent to the total star formation in all galaxies (Soifer & Nuegebauer 1998). 
Therefore the absence of them in UV selected samples does not introduce significant bias 
in the estimate of total starformation rate. But in the earlier universe of z > 1, this bias 
may be more significant. According to Le Floc'h et al. (2005), about more than 10% of star 
formation at z ~ 1 is due to ULIRGs. 
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